Design of mid-IR and THz quantum cascade laser cavities with complete TM photonic bandgap 
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We present the design of mid-infrared and THz quantum cascade laser cavities formed from planar photonic 
^ — ■ crystals with a complete in-plane photonic bandgap. The design is based on a honeycomb lattice, and achieves 

a full in-plane photonic gap for transverse-magnetic polarized light while preserving a connected pattern for 
f"^ , efficient electrical injection. Candidate defects modes for lasing are identified. This lattice is then used as a 

■ model system to demonstrate a novel effect: under certain conditions - that are typically satisfied in the THz 

range - a complete photonic gap can be obtained by the sole patterning of the top metal contact. This possibility 
greatly reduces the required fabrication complexity and avoids potential damage of the semiconductor active 
region. 

PACS numbers: 
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O |. Quantum cascade (QC) lasers are semiconductor laser sources based on intersubband (ISB) transitions in multiple quantum 
well systems^. Their emission wavelength can be tuned across the mid-infrared 2 ^ (mid-IR, 5 /jm < X < 24 /jm) and THzd (65 
^ \ /jm < X < 200 /jm) ranges of the electromagnetic spectrum. In the mid-IR, where they are becoming the semiconductor source 
O ■ of choice thanks to their good performance, the potential applications are chemical sensing, spectroscopy and free-space optical 
' ^ ! communications. In the THz range the most promising application is the imaging of concealed objects. Biological materials, 
£>v semiconductor chip packaging and clothing are all THz-transparent (while they are opaque at shorter wavelengths), making THz 
i-G ' imaging useful for security and medical applications. 

Most of the activity in the QC laser field has concentrated on edge-emitting lasers due to the intrinsic transverse-magnetic (TM) 
polarization of ISB transitions, and corresponding difficulty in implementing vertical-cavity surface emitting lasers (VCSELs). 
t— I \ Surface emission has been obtained, rather, by integrating second-order gratings on edge emitting devices^ or by replacing the 
J> ■ standard Fabry-Perot cavity with a photonic crystal (PC) resonator^. The latter solution, the application of photonic crystal 
technology to QC lasers*^, is particularly appealing because of the flexibility that it allows the designer. Two-dimensional (2D) 
photonic crystals can be used to create localized microcavity laser sources that can be built as two-dimensional arrays on a single 
| ' chip 11 ' 12 or for large-area, high-power, single-mode surface emitting laser sources^. 

Of particular importance in determining the properties of any photonic crystal structure is the effective refractive index contrast 
£ — ' attainable. The higher the index contrast, the stronger the optical dispersion of the photonic bands and the greater the ability to 
C"^) , localize, diffract, and reflect light within the photonic lattice. A high index contrast is therefore crucial for device miniaturization. 
""^2 ■ In planar waveguide devices, the presence of a full optical band gap in two-dimensions is beneficial (although not completely 
O \ necessary) in forming ultra-low-volume laser cavities with minimal optical loss. As is discussed in RefJ^, for 2D photonic 
" 55 ' crystals, TM optical band gaps are favored in a lattice of isolated high-e regions. Unfortunately, this configuration is incompatible 
£>V with an electrical injection device due to its non-connected nature and an alternate approach is required. 
fG In this paper we study the use of a connected honeycomb lattice for creating 2D photonic crystal QC laser structures. We begin 

G-r in Section IK] with a review of the properties of this 2D lattice via a planewave expansion (PWE) analysis. We show that a full 
2D optical bandgap for TM polarization can be obtained in this connected lattice, and study the localized resonant modes that 
form around a simple point defect in the lattice. In Sections|lII]and|lV]full three-dimensional (3D) finite-difference time-domain 
(FDTD) simulations are used to analyze the properties of the honeycomb lattice in two different vertical waveguide structures. 
Section |III]deals with mid-IR PC QC lasers. The real active-region and waveguide structure of a mid-IR surface-plasmon (SP) 
QC laser is considered. The high-index contrast in this waveguide structure is obtained by air holes that penetrate deeply into 
the semiconductor layers. Section [IV] focuses on THz PC QC laser waveguide structures. In the metal-metal waveguides of 
THz lasers we find that the effective index contrast, and therefore the photonic bandstructure, is strongly dependent on the 
waveguide thickness^. In particular, below a certain critical waveguide thickness, a full photonic gap can be induced by the sole 
patterning of the metal layers. This novel effect could be useful for the development of a variety of THz lasers, including PC 
surface-emitting lasers^, due to the simple fabrication requirements. 
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H. TWO-DIMENSIONAL ANALYSIS: PLANE WAVE EXPANSION 
A. The planewave expansion method 

The calculations of this section are based on a planewave expansion (PWE) method. It is a frequency-domain method that 
allows extraction of the Bloch fields and frequencies by direct diagonalization. The calculations were performed using the MIT 
Photonic Bands (MPB) simulation tool^. With the PWE method it is possible to compute the bandstructures and electromagnetic 
modes of perfectly periodic dielectric structures. Calculation of defect modes will be tackled in section HTCl using a super-cell 
approach. 

B. Photonic crystal structure and TM gap 

The 2D photonic crystal considered in this work is that of a honeycomb lattice, consisting of a hexagonal Bravais lattice with 
a basis of two air holes per unit cell (see Fig. [TJ a ))- The hole radius is r and the lattice constant is a. The reciprocal lattice is 
also hexagonal, with F, X and / the high symmetry points of the lattice (see Fig. QIb))- In the 2D simulations of this section the 
relative dielectric constant of the dielectric background is taken as e = 1 1 .22 (n = 3.35), corresponding to the index of refraction 
of most III-V semiconductor materials at frequencies below the energy gap. 




FIG. 1: Structure of the honeycomb lattice . (a) Two-dimensional honeycomb lattice with two "dielectric atoms" (A and B) per unit cell, 
ai = (a,0) and — ( a /2, v3a/2) are the principal lattice vectors (a = |ai | = |ao |), r is the hole radius, (b) Two-dimensional reciprocal space 
for the honeycomb lattice. The reduced Brillouin zone is shaded in blue. Reciprocal lattice vectors are superpositions of Gi = (0,4jl/V3<z) 
and G2 = (2n/a,2n/V3a). 




FIG. 2: Photonic band calculations for a 2D honeycomb lattice of air holes in a dielectric lattice (£ = 11.22). (a) Photonic bandstructure for 
TM polarization with r/a = 0.234. (b) Complete band gap (shaded blue region) as a function of r/a. A full TM gap exists only for r/a > 0.18 
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Fig. |2 a) shows the photonic bandstructure for a lattice with r/a = 0.234. For this r/a ratio a full TM gap is present, centered 
around a/A, ~ 0.225, with a width approximately 12% of the central gap frequency. The existence of a gap, however, depends 
strongly on r/a. In Fig|2b) we plot the gap-map for this lattice, which clearly shows that only above a critical value of circular 
hole size (r/a > 0.18) can a TM photonic gap be obtained. This is a common behavior of lattices of air holes; the larger the 
holes, the closer the lattice is to a system of isolated high-e regions, and the larger the TM photonic bandgap 9 . The advantage of 
the honeycomb lattice is that the TM gap opens when the hole size and porosity of the lattice are still reasonably small, allowing 
for more efficient electrical injection. 



C. 2D defect design: supercell method 



In this section we study, in 2D, the localized resonances that form around a point-like defect of the honeycomb lattice. The 
defect that we consider is obtained by removing a hexagon of six holes, as shown in Fig. [2 a). The defect modes are calculated 
with a supercell method using the same PWE solver. The presence of the defect cavity breaks the periodicity of the lattice, 
requiring the creation of a supercell over which the structure is assumed periodic. This supercell contains the defect cavity and 
is tiled in space. The Wigner-Seitz cell of the reciprocal lattice of the supercell will be correspondingly smaller than that of the 
underlying photonic lattice, resulting in a folding of the honeycomb lattice photonic bands (see Fig. Ob)). The localized defect 
states of the supercell appear as bands with almost no dispersion, lying inside the photonic bandgap of the honeycomb lattice. 
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FIG. 3: (a) Illustration of the 2D honeycomb lattice defect cavity supercell (background dielectric material shown as grey, air holes shown as 
white). The principle lattice vectors are ai, a2, where |aj | = [a^l = a. The defect region consists of removal of the central hexagon of air holes 
(central hexagon shown as a dashed line), (b) Folded TM photonic bandstructure for a 20-period supercell of the defect cavity shown in (a) 
with r/a = 0.24. The blue-lines correspond to the folded bandstructure of the honeycomb lattice. The red lines are the defect frequency levels. 
Localized defect modes: (c) hexapole and (d-e) dipole-like modes. 

A PWE calculated bandstructure of the supercell defect structure with r/a = 0.24 is shown in Fig. &b). In these calculations, 
a supercell consisting of 20 periods of the honeycomb lattice was used in order to obtain well localized resonant modes. Three 
defect modes lie within the complete band gap: two degenerate dipole-like modes and a hexapole-like mode (Figs. |3J c - e ))- 
These are the defect modes that are predicted for the honeycomb lattice through simple symmetry arguments^. Of the defect 
modes present in the bandgap, the dipole-like modes are more concentrated in the center of the defect region, thus giving them 
the most overlap with the QC gain material. As such, these are the cavity modes we will focus on in the 3D simulations described 
below. 

Our choice of hole radius (r/a = 0.24) in the above cavity design was made based upon a trade-off between the extent of the 
photonic bandgap and cavity mode localization, with that of the electrical and thermal resistance incurred in a semiconductor 
realization of such a structure. As already mentioned, the reduced connectivity of a photonic lattice makes electrical injection of 
a defect cavity more challenging. Electrical current in such structures is typically injected from the edge of the photonic crystal 
vertically through the device active region. Injection into the defect region of the cavity is a result of lateral current spreading 
into the center of the photonic crystal. The increased lateral resistance of highly porous photonic lattices significantly reduces 
the injection efficiency of the laser, resulting in added heating and reduced gain. Other injection geometries are possible, such 
as the use of surface-plasmon laser cavities in which the top metal contact extends over the entire photonic crystal and provides 
not only vertical waveguiding but also electrical injection. This geometry was used in the demonstration of the first electrically 
injected 2D photonic crystal microcavity laser—. Nonetheless, high-aspect ratio 2D photonic crystal lattices of high porosity 
pose a greater challenge to fabricate due to the reduced critical dimension size. It should be noted that further optimization of 
the honeycomb lattice can be obtained through use of air holes of a modified geometry. Lattices formed from air holes with 
a truncated circular cross-section were found to have larger photonic bandgaps for a given critical dimension than that of the 
standard lattice. In what follows, however, we focus on the honeycomb lattice with circular air holes. 
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HI. 3D-FDTD ANALYSIS OF MID-IR DEVICES 

The PWE method allows the rapid solution of structure eigenmodes in the frequency domain. However, for the 2D analysis de- 
scribed above one can only approximate the laser active region with an effective index of refraction. The objective of this section 
is to verify the design developed with the 2D model in Section[II] and to apply it to a realistic mid-IR QC laser structure. An ac- 
curate representation of the structure waveguide in the epitaxial direction will be taken into account within a 3D finite-difference 
time-domain (FDTD) approach. It is well known that the material system of choice for mid-IR QC lasers is InGaAs/AlInAs 
lattice matched to InR^. In addition, waveguides based on surface-plasmonsi 7 - have been shown to be advantageous for PC QC 
lasers in the mid-IRi^. This is therefore the model system that we will use for the 3D simulations. 

A. Mid-IR surface-plasmon waveguides for QC lasers 

Semiconductor diode lasers and conventional QC lasers rely on optical waveguides where a higher-index-core is sandwiched 
between thick cladding layers of lower refractive index, thus confining the light inside the active region stack. The small 
refractive index difference between the active core and the waveguide claddings results in a standard mid-IR QC laser which is 
typically 6 to 9 fjm thick. 




Distance (urn) Distance (,um) 

FIG. 4: Normalized intensity profile of the fundamental optical mode computed for a quantum cascade laser dielectric waveguide (a) and a 
surface-plasmon waveguide (b) designed for a wavelength of 8 /am. The shaded green areas indicate the stack of active regions and injectors. 
The vertical dashed line indicates the position of the device surface. In the dielectric case the thickness of the epitaxially grown material is ~6 
/jm. The origin of the abscissa is at the air-semiconductor interface. In the surface-plasmon case the thickness is instead ~ 3 /jm. The optical 
confinement factor is indicated by F. 

Maxwell's equations, however, also allow for another type of optical waveguiding based upon surface-plasmon&i- TM 
polarized electromagnetic guided modes exist at the interface of two dielectrics with opposite sign of the real part of their 
dielectric constants. Negative dielectric constants are typical of metals below the plasma frequency. Thus, guided surface- 
plasmon modes at a metal-semiconductor interface are a useful waveguiding solution for QC lasers^. This is due to the intrinsic 
TM polarization, i.e. normal to the layers, of intersubband transitions. SP waveguides need a smaller thickness of grown 
material, while yielding even larger optical confinement factors F (see Fig. |4j. 

The surface-plasmon damping along the propagation direction can be approximated with the following formula^: 

47t «,„•«/ 

where k m (n m ) is the imaginary (real) part of the metal index of refraction, is the real part of the semiconductor index of 
refraction, and A, is the wavelength. The 1 /k m X dependence of the propagation losses shows that SP waveguides are especially 
appealing at long wavelengths. In particular the l/k m factor, which pushes the field out of the metal region and reduces the 
losses, becomes very small at long wavelengths. A simple Drude model shows that the imaginary part of the index of refraction, 
for a generic metal, increases dramatically when moving from short (A. = 1-3 /jm) to long (A = 100-200 /jm) wavelengths. Recent 
advances have also shown that low-loss metallic waveguides can be implemented at mid-IR wavelengths 17 . We therefore employ 
this waveguide structure as a model system for the 3D numerical simulations of our QC laser structures. The corresponding layer 
sequence is displayed in Table|U together with the corresponding indeces of refraction that are used for the simulations. 
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TABLE I: Layer structure for the mid-IR SP QC laser 1 - modeled in the 3D FDTD simulations. Nominal operating wavelength is A, ~ 8 fjm. 



Material 



Doping (cm ) 



Thickness (/jm) 



Function 



air 

Au 
InGaAs 
InGaAs-AlInAs MQW 
structure 
InGaAs 

InP 



n.a. 

n.a. 
10 17 /10 18 
1.7 x 10 16 

5 x 10 16 

10 17 



4 

0.3 
0.05 
2.6 

0.5 

3 



1 

perfect conductor 
3.3/3.47 
3.374 

3.475 

3.077 



top cladding 
metal contact 
contact layers 
active region 

Buffer/lower cladding 
layer 
substrate 



While the vertical confinement is provided by the SP waveguide, the in-plane optical confinement is induced by the photonic 
crystal. A high index contrast can be obtained by penetration of the air holes through the top metal layer and deep into the semi- 
conductor waveguide structure. Intuitively, reduced scattering losses are obtained for a hole depth which overlaps a significant 
fraction of the guided mode energy 22 . In Fig. |3a) we show a schematic of the vertical cross-section of the photonic crystal 
SP waveguide structure, and in Fig. |3Jb) we plot the vertical mode profile of a localized defect mode (see Fig. |6|a)) of the 
honeycomb lattice of air holes with a hole depth of 4.7 /jm in the semiconductor heterostructure. At an operating wavelength 
of 8 /jm this results in a vertical mode overlap (energy overlap) with the air holes of almost 90%. It is clear from this example 
why surface-plasmon QC lasers are ideally suited to PC technology; their reduced thickness allows for a significantly shallower 
etch of the semiconductor material in comparison to conventional laser waveguide structures in which an etch depth of 10 /jm or 
more would be required in the mid-IR. 
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FIG. 5: (a) Surface-plasmon optical mode profile superimposed on a sketch of the vertical section of a mid-IR photonic crystal QC laser 
with air holes extending 4.7 jjm deep into the semiconductor, (b) 3D-FDTD calculated electric field intensity cross-section of a defect mode 
(top- view shown in Fig. [(Ja)) of the deeply patterned honeycomb photonic lattice. The vertical mode profile shows a similar behavior as the 
ID surface-plasmon simulation of (a), whereas the confinement in the in-plane direction is provided by distributed Bragg reflection of the 
photonic lattice. 



B. Defect cavity design 

The cavity investigated in this section is the same as in Section III CI It is obtained by removing a full hexagon of air holes 
from the PC lattice. As depicted in Fig. [5] we "etch" 4.7-/jm deep holes into the semiconductor laser structure, with the r/a ratio 
of the air holes set to 0.21. The properties of the different laser structure layers used in the simulation are given in TableQ] Full 
3D-FDTD simulation s 1 M 2 ' 22 were performed with an effective grid resolution of 100 nm, corresponding to roughly 20 points 
per wavelength in the high-index semiconductor material at an operating wavelength of 8 /jm. The use of perfectly conducting 
boundaries for the top metal contact layer is a small approximation at 8 /jm due to the much higher plasma frequency of the Au 
metal contacts used in practice. 

A simulation of a photonic crystal cavity with n x = l periods in the i-direction and n y = 4 periods in the y-direction yields 
the three localized defect modes shown in Fig. HJa-c). These modes are the same modes as found in the 2D analysis of Fig. [3] 
Comparison of the 2D and 3D mode properties are summarized in Table [II] The normalized frequencies of the defect modes 
lie within the TM-like bandgap of the 3D honeycomb lattice waveguide structure, and are centered around a/X ~ 0.2. At a 
wavelength of 8 /jm this corresponds to a lattice constant of approximately a = 1.6 /jm, a hole radius of r = 0.34 /jm, and a 
minimum dimension given by the gap between nearest neig hbor holes of a (I /y/3-2r/a) = 0.25 /urn. 

Effective quality factors associated with in-plane (Qn), topside vertical (Q t ), and bottomside vertical (Qt,) radiation losses are 
also calculated for the 3D FDTD simulations. The topside vertical effective Q-factor for all of the defect modes is estimated at 
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TABLE II: Comparison between 2D (PWE) and 3D (FDTD) simulated defect modes. 



Parameter 



x-dipole mode 



y-dipole mode 



hexapole mode 



PWE a/X 
FDTD a/X 

Q\\ 

Qi 
Qb 



0.2115 
0.196 
326 
1.2 x 10 7 
92 



0.2115 
0.197 
288 
9x 10 6 
93 



0.2162 
0.205 
169 
1.1 x 10 7 
114 



> 10 , due to the fact that the cavity modes lie predominantly below the light cone of the top air-cladding . Qb on the other 
hand, can be increased only slightly from that reported in Table HI1 through deeper etching of the air holes, limited mainly by 
the high-index of the semiconductor substrate and the localized nature of the cavity. Due to the TM-like polarization in-plane 
bandgap provided by the honeycomb lattice, an increase in the number of PC mirror periods to a value above 10 effectively 
eliminates in-plane radiation loss, and the Q-factor of the cavity modes is limited by both radiation into the substrate as well as 
material absorption loss due to ohmic heating in the metal surface-plasmon layer (for wavelengths in the mid-IR, and a Au metal 
surface-plasmon layer, material absorption limits the Q-factor to approximately 10 3 — ). 
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FIG. 6: In-plane mode profiles (£-) for the (a) x-dipole, (b) j-dipole, and (c) hexapole defect modes of the surface-plasmon vertical waveguide 
structure with an "etched" hexagonal defect cavity in the honeycomb lattice. 



IV. ANALYSIS OF THZ DEVICES 



Our FDTD analysis of mid-IR PC QC lasers examined structures in which the honeycomb lattice was patterned through the 
QC semiconductor active region. This is necessary in order to induce a strong index contrast. However, a different solution 
is possible if the optical mode is highly confined in the vertical direction, as is the case for THz QC lasers with metal-metal 
waveguide structure o 24 ' 25 . In this case, and under certain constraints related to the structure thickness, the sole patterning of the 
top metal contact is capable of inducing a complete photonic bandgap. These metal-insulator-metal (MIM) structures 2 ^, as we 
will refer to them here (see Figs. |7]and[8]l, will be the focus of our time domain modeling in the following sections. The appeal 
of this design lies in the fact that only the top metal layer needs to be patterne d 27,28 , which both simplifies the fabrication process 
and improves the efficiency of carrier diffusion in the active region. 



A. MIM structures: waveguides for THz QC lasers 



Metal-metal waveguides^ are usually employed for QC lasers in the THz rang o • 5 because they can provide almost unity 
optical confinement factors, and simultaneously, relatively low waveguide losses. The active laser core is sandwiched between 
two metal layers (typically Ti/Au or Ge/Au/Ni/Au) which act as surface-plasmon layers (Fig. Ha)). The two surface-plasmon 
modes of both the top and bottom metal-core interfaces become coupled and form two guided modes, one of even parity and 
one of odd parity (assuming a nearly symmetric metal-core-metal structure, and where the vector parity is determined by the 
symmetry of the magnetic field). The dispersion diagram of a double-metal waveguide with active layer core thickness of L a = 1 
fjm is displayed in Fig. Ui e )- F° r sub-wavelength core thicknesses the odd parity surface mode is cut-off, while the even parity 
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mode exists all the way down to zero frequency. The electric and magnetic field profiles of the even parity surface mode at a 
wavelength of A,q = 100 /jm are shown in Fig. [7Jb-d), showing the near-unity confinement factor of the fields in the core region. 




FIG. 7: (a) Schematic showing the cross-section of a MIM waveguide (metal layers are assumed semi-infinite), (b) E z normal electric field 
component, (c) E x longitudinal field component, and (c) H v magnetic field plot of the even parity guided surface mode at X = 100 /am for an 
active region core thickness of L a = 1 /jm. (e) Dispersion diagram of the double-metal (Au) waveguide structure (L n = 1 /jm). The core is 
modeled with a constant refractive index of n a = 3.59, while the metal layers are modeled using a Drude-Lorentz model for Au (background 
dielectric constant Z\, = 9.54, plasmon frequency a p = 1.35 x 10 16 rad/s, and relaxation time 8 x 10~ 15 s). 



As has been noted by other authors^, a reduction in the dielectric core region thickness of MIM structures results in an 
increase in the in-plane wavevector of the even parity guided mode for a particular frequency. This is a result of the "pushing" 
of the electromagnetic energy into the metal cladding regions with decreased core thickness. Of import to the current work is 
the ability to create a large effective index contrast through patterning only of the metal layers. In such a patterned structure the 
guided mode sees two effective indeces as it propagates, the large effective index of the guided surface-plasmon mode where 
the metal is left intact, and a lower effective index, more delocalized mode of the dielectric core and air-cladding where the 
metal is removed. That this picture is in fact valid, and can be used to great effect, is demonstrated below where we analyze the 
bandstructure and localized resonances of double-metal waveguides with only a surface-patterning in the top metal contact. 



B. Patterned MIM waveguide bandstructure analysis 



6°0 0°0 

o o o o o o 

o°o° rh?o°o 

°o9o o°o° 

o%9o°o°o°o 
o o o o o o 



fb) 


Air 


perfect conductor 




QC active region 

(n 3 =3.59) 


perfect conductor 



t 5ir = 25nm 



t m = 1 urn 



fm = 1 |im 



FIG. 8: Schematic of the honeycomb lattice and MIM structure used in FDTD modeling, (a) Top view of the lattice with a central defect 
consisting of the removing of the central hexagon of air holes, (b) Cross-section of the simulated structure. Only the top metal layer (perfect 
conductor) is patterned. The displayed layer thicknesses are for an operating wavelength of Xq = 100 and a normalized frequency of 
a/Xo = 0.17 within the bandgap of the honeycomb lattice. 



The model MIM structure that we will consider here, shown schematically in Fig. Ob), consists of the following sequence 
of layers (from bottom to top): a bottom unpatterned metal contact, a QC active region dielectric core (modeled with a uniform 



8 



index of n a = 3.59 for simplicity^ -), a top patterned metal contact, and a region of air above the structure. The metal layers are 
assumed to be perfect conductors; loss due to absorption in the metal regions is discussed in section HV D| Mur absorbing bound- 
ary conditions^ are used to model radiation loss out of the structure. For each of the THz FDTD simulations described below, a 
grid resolution of 58 points per lattice constant a of the honeycomb lattice was used. At a nominal operating wavelength of 100 
jum, and for normalized frequencies a/Xo ~ 0.17 within the bandgap of the honeycomb lattice (see below), the lattice constant 
is a = 17 /jm. At this wavelength and normalized frequency (hereinafter the nominal operating conditions), the corresponding 
spatial resolution is 0.3 /jm (roughly 95 points per wavelength in the dielectric core material), and the MIM layer thicknesses are 
as shown in Fig. |8jb). 




□ X J □ 

FIG. 9: FDTD calculated bandstructure of a patterned MIM waveguide with a normalized dielectric core thickness of L a = L a /a = 0. 176 and a 
normalized air hole radius of r/a = 0.25 in the top metal contact. Metal layers are simulated as perfect conducting boundaries. At an operating 
wavelength of Xq = 100 /mi within the photonic bandgap (ci/Xq = 0.17), the corresponding physical sizes of the patterned double-metal 
waveguide are {a, r, L a } = {17, 4.25, 3}pm. 

The band diagram for surface-patterned MIM waveguide structures with various dielectric core thicknesses were computed 
to confirm the existence of a complete in-plane bandgap, and to examine the dependence of the bandgap frequency width 
on structure thickness. The simulation volume consists of a single unit cell of the honeycomb lattice with Bloch periodic 
boundary conditions applied in the plane of periodicity and Mur boundary conditions in the top and bottom directions (normal 
to the metal and semiconductor layers). Figure [9] shows the band diagram for a structure with a normalized core thickness of 
L a = L a I a = 0.176 and a normalized air hole radius of r/a = 0.25 in the top metal layer. The corresponding physical core 
thickness is L a = 3 under nominal operating conditions (Xq = 100 /jm, a/Xo = 0.17). In the band diagram, truly guided 
modes lie below the air light-line, while leaky modes lie above it (in this case light can only leak into the top vertical direction 
where patterning of the metal has been applied). The bottom-most band is a surface wave, lying just below the light-line of 
the active region core. It is akin to the even parity guided surface-plasmon mode of the symmetric (unpatterned) double-metal 
waveguide studied above. The higher-lying frequency bands constitute zone-folded versions of this surface wave. An in-plane 
photonic bandgap exists for guided modes of this structure between normalized frequencies a/Xo = 0.165 — 0.18. Due to the 
extreme thinness of the dielectric core, the higher-order vertical modes of the MIM waveguide are not present in this diagram, 
but lie at much higher frequencies (a/Xo > 0.75). As a result, the photonic bandgap is a true full in-plane bandgap for guided 
resonances of the MIM structure, not just for a single polarization or mode symmetry. Note, the modes that lie along the air 
light-line in the diagram are radiation modes, predominantly localized in the air regions above the MIM structure. 

In Fig. [10] we plot the dispersion of the bands defining the (lowest) photonic bandgap between the high symmetry X and J 
points of the honeycomb lattice for varying dielectric core thicknesses. The bandgap is seen to shrink with increasing dielectric 
core thickness, closing for a normalized active region thickness greater than L a = 0.294 (nominal physical thickness L a =5 /jm). 
Mode field plots (see Fig. ITTb at the X and J points indicate that the high and low frequency modes defining the bandgap are 
of mixed type: the low-frequency "dielectric" or "valence" band mode is predominantly a surface wave at the interface of the 
top patterned metal surface, whereas the high-frequency "air" or "conduction" band mode is predominantly a surface wave of 
the bottom unpatterned metal interface. Additionally, the high-frequency gap mode sits largely beneath the unpatterned regions 
of the top metal boundary and the low-frequency mode resides below the patterned areas with more energy residing in the air 
cladding. The photonic bandgap in this case is a result of differences in these two surface waves. This simple picture, although 
correct, betrays its complexity for two reasons: (i) perfect metal conductor boundaries were used in the simulations, which do 
not support surface waves when continuous and flat, and (ii) the patterned air holes are much smaller than the wavelength in air 
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FIG. 10: (a) Valence and conduction band dispersion between the X and J points for varying active region core thickness (L a = 
{0.088(d), 0.176(x),0.294(o), 0.412(D)}; nominal physical thicknesses L a = {1.5(d), 3(x),5( >),7(D)} /mi). The symbols correspond to 
data points from FDTD simulations and the lines are guides to the eye. The normalized hole radius for all simulations was fixed at r/a = 0.25. 
Note that the complete band gap shrinks as device thickness increases, and is closed for L a > 0.294 (L a > 5 /jm). 

of the guided mode. That a surface wave can indeed exist at the interface of a patterned perfect conducting layer was pointed 
out in Refi^, where it was shown that it is precisely the cut-off frequency of the air holes in the perfect conductor which sets 
the effective plasma frequency of the layer. The thinner the dielectric core layer, the larger the fraction of energy that resides 
at the top patterned metal interface, and a greater degree of mode derealization that occurs into the air-cladding where the 
top contact is removed, both effects which increase the effective index contrast of the two surface waves and thus the photonic 
bandgap. It should be noted that the use in our model of perfect conducting boundaries as opposed to a Au metal contact layer, 
for instance, does not affect the photonic bandstructure shown in Fig. [9] and Fig. [10] as the plasma frequency of Au is much 
greater than the THz frequencies under consideration. Further discussion and analysis of photonic bandgap structures formed 
from surface patterned metal and perfectly conducting layers will be presented in a future study^. In practice, laser structures 
with thicknesses much below that of several microns may be difficult to fabricate and incur insurmountable optical losses (see 
below for discussion), but the trend shown in Fig. [10] clearly indicates that such structures would support complete bandgaps 
that are quite substantial, approaching 15% of the gap center frequency for L a = 1.5 /jm at an operating wavelength of A,o = 100 
/jm. 



Ultimately, one would like to utilize the strong photonic bandgap effects present in the patterned MIM waveguide structures 
studied above to, for instance, form localized microcavity laser resonators. In order to demonstrate the effectiveness of the top 
metal patterning in this regard, we have also simulated the localized resonances that are formed within and around point defects 
of the honeycomb metal pattern studied in the previous sub-section. In Fig. [TZl a) we show the E z field plot of the localized 
y-dipole-like^ mode that forms around a point-defect consisting of the removal of the central 6 holes in the top metal (see Fig. 
[HJa)) of a MIM waveguide with an active region thickness of L a = 3 /jm (L a = 0.176). The corresponding (in-plane) spatial 
Fourier transform of the mode is plotted in Fig. [T2l b). showing that the mode is localized at the kx points in reciprocal space, 
and consistent with a "donor" mode formed from the "conduction" bandedge at the X-pointii. This mode is identified as the 
same dipole-like mode of the mid-IR surface-plasmon simulations of Fig. |6jb) in section!]]]] where in those structures a deep air 
hole patterning of the semiconductor dielectric material was used to create a large index-contrast. 

For the fixed waveguide core thickness of 3 /jm the y-dipole-like mode of Fig. [T2ja) is seen to be highly localized in real- 
space, as expected from the bandgap simulations of Fig. [10] Additional confirmation of the in-plane localization of the defect 
mode due to an in-plane photonic bandgap effect can be obtained by studying the scaling of the in-plane radiation losses as a 
function of the number of periods of the PC patterning surrounding the central defect. The total quality factor due to radiation 
loss is determined by 



C. Characterization of defect modes 
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FIG. 1 1 : High-frequency (a-c) and low-frequency (d-f) gap mode field plots (E z ) at the X -point of the honeycomb lattice for the patterned 
double-metal waveguide structure of Fig. [lO]with L a = 0.294 {L a = 5 /mi). The in-plane field plots of (a,d) correspond to a section through 
the core of the MIM structure just below the patterned top metal contact. The dashed white lines in (a,d) and (b,e) indicate the position of the 
sections used in the field plots of (b,e) and (c,f), respectively. 

where (2|| and (2± quantify the in-plane and out-of-plane radiation losses^. The out-of-plane quality factor is found to be > 10 5 
for these structures; thus the dominant radiation loss occurs in the in-plane direction. Fig. [T3] shows a plot of the g-factor as 
a function of the number of surrounding PC periods. The radiation Q-factor increases exponentially with period number as 
expected for a bandgap mode, with a cavity consisting of only 12 x 6.5 periods able to sustain a mode with radiation Q > 10 3 . 



D. Effective mode volume and metal absorption losses 



We have shown above that in double metal structures the frequency-width of the photonic bandgap of an honeycomb photonic 
lattice increases when the structure is thinner. It implies, intuitively, that the photonic crystal becomes increasingly effective in 
confining light when the waveguide core thickness is reduced. To quantify the modal localization versus core region thickness 
we calculated the effective mode volumes V e ff of the y-dipole-like mode of Fig. [12] for various active region thicknesses (L a ). 
Veff is defined as: 



Je|E| 2 ^ 
max [e | E | 2 ] 



(3) 



where max[..] denotes the maximum value of the argument. The results are reported in Table Hill The decrease of V e s with the 
active core region thickness L a is super-linear, as can be seen in the third column of Table HIT1 in which V e g/L a is tabulated. The 
super-linear decrease in V e ff with active region thinning is a clear indication that the photonic crystal formed from patterning of 
the top metal surface provides more effective in-plane confinement for thinner waveguide cores. 

FDTD simulations to this point have approximated the metal contacts defining the MIM waveguide structures with perfect 
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FIG. 12: In-plane mode profile of the y-polarized dipole-like defect mode for a point defect consisting of the removal of the central 6 air-hole 
patterns in the top metal contact (see Fig. [8}. (a) Real space profile of E- and (b) magnitude of the spatial Fourier transform of E z for an 
in-plane cut through the middle of the MIM waveguide. 
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FIG. 13: Increase of Qu with the number of PC periods for an MIM structure with L a = 3 /jm (L a = 0.176) thick active region. The diamonds 
show the results of FDTD simulations and the solid line is a least squares fit to these results. The PC cavities were chosen to be approximately 
square, so that the number of lattice periods (n x x n y ) used in the FDTD simulations were 3x2, 5x3, 7x4, 9x5 and 11x6 (referring to Fig. 
[Ha), a "period" in the x direction is taken as a, whereas in the v direction it is \/3a). The abscissa of the plot refers to the number of periods 
(n x ) in the x direction. 



conductor boundary conditions, neglecting ohmic heating loss in the metal layers. Additionally, the optical losses within con- 
ventional THz QCLs include a contribution from free-carrier absorption, the level of which can be minimized through reduction 
of the thickness and doping levels of the semiconductor contact layer a 24,25 ' 30 . Given that the thin MIM waveguide structures 
studied here have a large fraction of mode energy at the metal surfaces, we will concern ourselves here with estimating the opti- 
cal losses in the metal layers, with free-carrier absorption and other losses contributing at a level similar to that in conventional 
THz QCL structures. Incorporation of "real" metals into the 3D FDTD simulations would require a prohibitively large amount 
of memory and/or simulation time due to the very short skin depth at THz frequencies. Fortunately, an estimate of the ohmic 
heating losses in the metal can be obtained for the patterned MIM waveguide structure by considering the analytically tractable 
unpatterned MIM structure. 

Specifically, we approximate the losses due to absorption in the metal contacts through an analysis of the fundamental, 
even parity, coupled-surface-plasmon mode in a simple planar ID double-metal waveguide structure 2 ^. The complex dielectric 
constant of the Au metal layers in the THz (X = 100 /jm) are taken from RefM- (e Au = -1.06 x 10 5 + /1.93 x 10 5 ). The 
refractive index of the QC active region core (which is assumed to be homogeneous and lossless) is taken to be 3.59. The fields 
and complex propagation constant are determined through numerical solution of the dispersion relation. The imaginary part of 
the propagation constant, which determines the modal loss per unit length down the waveguide (see TableHIIb. is converted to a 
loss per unit time and related to an effective Q due to metal absorption for comparison with the radiation losses plotted in Fig. 
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FIG. 14: Effective g-factor due to metal waveguide loss (Q m , blue line) and effective energy velocity index (he = c/Vg, red line) versus active 
region thickness (L a ) for a ID MIM waveguide with Au metal guiding layers at Xq = 100 /jm. 
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2v E Im(|3) 

where Ve is the global energy velocity of the mode in the waveguide and coo the angular frequency of the mode. The energy 
velocity is computed from 



v E = t,' , , (5) 
f(W) dz 

where (5) is the time-averaged Poynting flux, (W) the time-averaged energy density valid for a lossy, dispersive medium—, and 
the quantities are integrated along the transverse width (z-direction) of the waveguide. Fig. [14] shows the variation of Q m and 
c/ve versus the width of the core region. This behavior is due to the penetration of the electric field into the lossy metallic layers 
when the active region is made thinne r 25 ' 29 . The metal-absorption-limited propagation loss per unit length (oc m ) and Q-factor 
(Qm) for active core region thicknesses of L a = 1.5, 3, 3.5, 5, and 7 /jm are tabulated in Table HIT1 along with the effective mode 
volume data. From this estimate we see that metal absorption loss will dominate the in-plane radiation losses of the y-dipole-like 
mode of Fig. Q~2]for defect cavities of 9 x 5 periods or more. 

In order to connect this work with that of more conventional THz QC laser structures and measurements, and to determine 
whether such thin double-metal waveguide structures are at all practical, one must relate the calculated Q-factors of Fig. [T3land 
Fig. [l4]to per unit length exponential loss coefficients. The typical gains that can be achieved in THz QC lasers at A, « 100 
jjm are ~ 20 — 40 cm" — . These values refer to modes whose energy velocity is approximately equal to the phase velocity of 
the bulk material. In thin metal waveguide structures the energy velocity of the plasmon mode can be much reduced relative to 
that for a mode in a thick waveguide structure (see the dispersion of Fig. lite)); however, for the MIM waveguide thicknesses 
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considered here (L a > 1.5 /jm) the energy and core phase velocity are nearly equal and the metal absorption loss coefficients, oc m , 
of Table HHI can be directly compared to measured gain coefficients of conventional laser structures. Such comparisons indicate 
that for all but the thinnest waveguide structure (L a = 1.5 fjm), laser action should be achievable. 



V. CONCLUSIONS 



In conclusion, we have described photonic crystal guided wave structures tailored for quantum cascade lasers, in which a full 
TM photonic bandgap is present and lateral electrical injection is possible. Laser cavity designs in both the mid-IR and THz 
ranges of the electromagnetic spectrum were considered. We identified candidate defect modes for lasing, and calculated their 
mode profiles and (2-factors. In addition, we have shown that it is possible - under appropriate conditions - to implement high- 
contrast photonic structures via the sole patterning of the top metal layer in THz QC lasers featuring double-metal waveguides. 
Such photonic crystal structures are of particular interest for laser structures due to the negligible damage they introduce into the 
active semiconductor layers and their reduced fabrication complexity. 
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